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ABSTRACT 

We highlight how the mass-radius relation of cluster-forming cores combined with an 
external tidal field can influence infant weight-loss and disruption likelihood of clusters at the 
end of their violent relaxation, namely, when their dynamical response to the expulsion of their 
residual star-forming gas is o ver. Specifically, building on the cluster A^-body model grid of 
iBaumgardt & Krou pa ( 2007), we study how the relation between the bound fraction of stars 
staying in clusters at the end of violent relaxation and the cluster-forming core mass is affected 
by the slope and normalization of the core mass-radius relation. Assuming mass-independent 
star formation efficiency and gas-expulsion time-scale Tge.xp/ Across and a given external tidal 
field, it is found that constant surface density cores and constant radius cores have the poten- 
tial to lead to the preferential removal of high- and low-mass clusters, respectively. In contrast, 
constant volume density cores result in mass-independent cluster infant weight-loss, as sug- 
gested by some observations. These trends result from how core volume density and core 
mass scale with each other. Infant weight-loss is quantified for cluster-forming cores with ei- 
ther number density n^.core — 6 x 10 era" 3 , or surface density T. core ~ 0.5 g.cm~ 2 , or radius 
fcore = 0.3 pc. Our modelling includes predictions about the evolution of high-mass cluster- 
forming cores (say m core > 10 s Mq), a regime not yet covered by the observations. We show 
how, for a given external tidal field, the core mass-radius diagram constitutes a straightfor- 
ward diagnostic tool to assess whether the tidal field influences the fate of clusters after gas 
expulsion. 

An overview of various issues directly affected by the nature of the core mass-radius 
relation is presented. In relation with the tidal field impact, these are the evolution of the 
cluster mass function at young ages (i.e. over the first ~ 30 Myr), and our ability to reconstruct 
the star formation history of galaxies from their cluster age distribution. Independently of 
the tidal field impact, the slope and/or normalization of the cluster-forming core mass-radius 
relation also influences the mass-metallicity relation of old globular clusters predicted by self- 
enrichment models, and the duration of cluster violent relaxation. 

Finally, we emphasize that observational mass-radius data-sets of dense gas regions must 
be handled with caution as they may be the imprint of the molecular tracer used to map them, 
rather than reflecting cluster formation conditions. 

Key words: stars: formation — galaxies: star clusters: general — galaxies: evolution — stars: 
kinematics and dynamics 



1 INTRODUCTION 

Modelling the early evolution of star cluster systems provides 
crucial insights into cluster formation physics through a compari- 
son between predicted and observed correlations and distribution 
functions of individual cluster properties. Over the past 30 years, 
considerable efforts have been dedicated to modelling cluster 
violent relaxation, i.e. cluster evolution after residual star forming 



gas expulsion (e.g. Tutukov 


1 19781: Iffillsl 19801: Mathieu 19831: 


Lada, Margulis & Dearborne 


Il984l: iKrouoa, Aarseth & Hurlev 1 


1200 ll; iGever & Burkert 


200l|; Goodwin & Bastianl |200fj|: 



l2009h. In 



partic ular, cluster iV-body model grids (e.g. Baum gardt & Kroupa I 



l2007h allow to model entire star cluster systems while browsing 
the parameter space extensively. Initial conditions of star cluster 
systems constitute crucial ingredients of their time-evolution 
modelling, and the mass-radius relation of cluster-forming cores 
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is therefore an issue at the forefront of the physics of both cluster 
formation and cluster evolution. 

In an influential study of Gala ctic molecular c louds and of the 
density enhancements they contain. lLarson I dl 98 ll) find that molec- 
ular clouds and their cores are in approximate virial equilibrium 
and have a near constant surface density. These properties were 
then included by lHarris & Pudritz ] d 19941) in their model of for- 
mation of old globular clusters, whose birth sites are identified as 
the dense cores of 'supergiant molecular clouds' at the early pro- 
togalactic epoch. Constant surface density cluster-forming cores 
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which contrasts 



show a strong mass-radius relation r core oc m core , 
with the absence of a clear mass -radius relation for gas-free star 
clusters, regardless of th e ir age llZepf et al.l 1 19991 : lLarsen 1 l2004 
IScheepmaker et al]|2007l) . iKroupa 1 1 20051) therefore suggests that 
cluster-forming cores themselves have uncorrelated masses and 
radii. Such an hypothesis has a bearing on the time-evolution of 
the mass function of clusters from the embedded phase to the 
end of violent relaxation. Owing to their deeper potential well, 
massive cores undergo more adiabatic g as expulsion and, there- 
fore, retain a greater fraction of their stars. IKroupa & Boilv I d2002h 
show that this can account for the for mation of a turnover in 
the cluster mass function a t young ages. iBaumgardt et al.l j2008t) 
and IParmentier et alj feOQS 1 ) further that hypothesis and show that 
cluster-forming cores with constant radii indeed produce features 
(flattening or turnover) in the cluster mass function provided 
that the star formation efficiency (SFE), assumed to be mass- 
independent, is not higher than 30 per cent. These studies mostly 
aim at explaining the prominent and universal tur nover characteriz- 
ing th e mass function of old globular clusters (see lAshman & Zepfl 
1 1998L and references therein). 

The mass function of young star clusters in the present-day 
Universe is reported to be a featureless power-law of spectral in - 
dex of about -2 (e.g. IZhang & Fall 1 1 19991 ; lLada & Ladall2003h . 
irrespective of the cluster age (say, 1, 10 or 100 Myr). This im- 
plies that cluster infant weight-loss (i.e. the gas-expulsion-driven 
cluster star-loss) is independent of the emb edded-cluster mass. Re- 
cently, iFall. Krumholz & Matznerl d2010l) have estimated that the 
SFE required for a cluster-forming core to expel its gas via stellar 
feedback is mass-independent if cluster-forming cores have a con- 
stant surface density. Their model assumes that the gas expulsion 
time-scale in units of a core crossing time, XcExpl Across, is constant. 
While preserving the shape of the cluster mass function at young 
ages, such a finding leaves unanswered the question of why young 

star clusters are deprived of a significant mass-radius relation if the 

1 12 

mass-radius relation of their progenitors scales as T core x Hlcore • 

In this contribution we are addi ng one more piece to 
this intriguing puzzle . Previous studies dKroupa & Boilv 1 1 20021 ; 
IParmentier etalj|2008l ; frail. Krumholz & Matznerll2010l) have ig- 
nored the influence that an external tidal field may exert upon clus- 
ters experiencing violent relaxation. Most cluster stars venturing 
beyond the cluster tidal radius become unbound field stars. Q There- 
fore, as a cluster expands following gas expulsion, its infant weight- 
loss and its likelihood of disruption are partly governed by how 
deeply the embedded cluster sits within its limiting tidal radius, that 
is, how severe tidal overflow due to cluster expansion is. lGoodwinl 
d 19971) performs A'-body simulations highlighting this effect: in- 
fant weight-loss of otherwise identical model clusters is stronger 



1 Stars on highly excentric orbits may experience transient passages be- 
yond the tidal radius and re-integrate into the cluster thereafter. However, 
such stars are expected to be rare. 



closer to the Galactic centre by virtue of the stronger tidal field and 
hence smaller cluster tidal radius (his fig. 3). IBaumgardt & Kroupa I 
d2007l) quantify this effect by the ratio of the half-mass radius r/, to 
the tidal radius r, of the embedded cluster. If r^/r, = 0.01, the clus- 
ter has much space in which to expand following gas expulsion, 
the tidal field impact is low and cluster infant weight-loss is solely 
driven by the SFE and the gas expulsion time-sca l e Tqe xd I ' Across- 
In contrast, if r h /r t = 0.20. lBaumgardt & Kroupa I d2007h find that 
protoclusters are mostly disrupted. We refer to rp,/r t as the tidal 
field impact. As we shall see in Section [2] not only is it related 
to the external tidal field, it also depends on the embedded clus- 
ter mass and size and, therefore, to the cluster-forming core mass- 
radius relation. It is therefore crucial to quantify to what extent the 
mass-radius relation of cluster-forming cores combined to an exter- 
nal tidal field influences the early evolution of star cluster systems. 
In this introductory paper, we focus our attention on the fraction 
of stars which remains bound to their parent clusters at the end of 
violent relaxation. 

The outline of the paper is as follows. Section [2] investigates 
how different cluster-forming core mass-radius relations (constant 
radius, constant volume density and constant surface density) influ- 
ence the cluster bound fraction as a function of core mass. We also 
show how the mass-radius diagram of cluster-forming cores can be 
used to estimate whether an external tidal field influences cluster vi- 
olent relaxation or not. To what extent mass-size data-sets of dense 
molecular gas regions can help us constrain the cluster-forming 
core mass-radius relation is the topic of Section [3] In Section [4] 
we comment about the importance of the core mass-radius relation 
regarding crucial issues such as cluster infant mortality /weight-loss 
as a function of cluster mass, and the reconstruction of the star for- 
mation history of galaxies based on their surviving clusters. We 
conclude in Section[5] 



2 CORE MASS-RADIUS RELATION AND EXTERNAL 
TIDAL FIELD 



At the end of its violent relaxation (age ~ 30 Myr; see Section |43t 
the mass of a star cluster is 



m c i =F bound .SFE.m c 



(1) 



with 

wicore mass of the cluster progenitor core and F DOlmd the 
fraction of stars remaining bound to the cluster at the end of violent 
relaxation. SFE is the 'local' star formation efficiency, namely, the 
mass fraction of gas turned into stars at the onset of gas expulsion. 

Cluster infant weight-loss, 1 — F ooulu j, is a sensitive func- 
tion of the SFE at the onset of gas expulsion and of the gas 
expulsion time-scale expressed in unit s of a clus te r-forming 
core c rossing-time, TgExpl Across ( e -g- Hills I 19801: Mathieu 



1983 1: lLada. Margulis & Dearbome I Il984l ; iGever & Burkert 
200 ll) . The higher the SFE, the slower the gas expul- 
sion time-scale XcExpl Across, the higher the bound fraction 

Fbound(SFYi,T(jE X p/T cross ). 

Formally, the bound fraction F aound depends on the 'ef- 
fective' star formation efficiency (eSFE , IVerschuerenl Il990l : 
iGoodwin & Bastian ||2006| ; iGoodwin ifioolj) . rather than on the lo- 



cal SFE. That is, F^ 



■ F bound {eSFE,T GExp /z cross ). The eSFE 



incorporates how far from virial equilibrium the cluster is at 
the onset of gas expulsion. If the stars and gas are in virial 
equilibrium at the onset of gas expulsion as a result of, for 
instance, a several-crossing-time time-span between star forma- 
tion and gas expulsion, the eSFE is simply the 'local' SFE and 
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FboundieSFE ,T GExp / tcross) = F bomd (SFE, %GExp I Across)- This is 
the approach we adopt in this contribution since the assump- 
tion of virial equilibr i um u nderpins the A'-body model grid of 
iBaumgardt & Kroupa I j2007l) (and the vast majority of other stud- 
ies dedicated to cluster gas expulsion). However, if the dynami- 
cal state of the newly-formed stars at the onset of gas expulsion 
is 'cold', then the eSFE is higher than the local SFE, and the 
bound fraction F 'bound becomes larger than shown in Figs. [3] [6] 
and [7] (and the opposite if stars are in a 'hot' dynamical state). 
Cluster models in which stars are n ot in virial equilibrium at gas 
expuls i on have been i nvestig ated by Lada. Margulis & Dearbornel 
< ll984l) . IVerschuerenl dl99dh and lGoodwinl l ll997h . Note however 
that if the cold collapse at gas expulsion stems from stars form- 
ing out of a contracting pre-cluster core, this implies that the star- 
formation activity in the pre-cluster cloud core would need to be 
synchronised to occur within a time shorter th an the core cro ssing- 
time and so SFE ~ eSFE is most likely (see iKroupa 1120081 for a 
discussion). 

In addition to the SFE and gas expulsion time-scale, the bound 
fraction F oomc i may also depend on the tidal field impact o,/r r , i.e.: 

Fbound = Fbound{SFE ,T CExp / 'T c rois,Oi / ''"') i (2) 

an effect mapped by IBaumgardt & Kroupa I J2007I) by means of 
cluster A'-body modelling. Stronger tidal field impacts r^jr, lower 
the bound fraction Fbound- To provide a clear understanding of 
how the cluster-forming core mass-radius relation affects the bound 
fraction through the tidal field impact, in what follows, each simula- 
tion is assigned a given local SFE, a given gas expulsion time-scale 
and a given external tidal field. That way, any Fb olmc t variation will 
necessarily result from varying the core mass-radius relation. 



2.1 Fiducial model: SFE=0.33, Zgexp — Across 

We adopt SFE = 0.33 dLada & Lada 1 120031) and x CExp ~ x cr0 s S 
dKrumholz & Matzner I2009I . see also Section HTt as the local SFE 
and gas expulsion time-scale of our fiducial model. The parameter 
space in terms of SFE and Xqexp I ' Across ls explored more widely 

later in this section. 

At this stage, we note that IBaumgardt & Kroupa (2007) 
model cluster gas expulsion as an exponential decrease with time 
of the cluster gas mass m gas (see their eq. 3): 



(t) 



,(0)« 



(3) 



Therefore, Tm, which they define as the gas-expulsion time-scale, 
is actually the e-folding time of the gas expulsion process. It 
corresponds to the time when the residual gas mass fraction is 

= 0.37 of its initial value. In all our models, we define the 
gas-expulsion time-scale %GExp rather as the time-scale over which 
the cluster expels the entirety of its residual gas. Pri or to using the 
A'-body model grid of IBaumgardt & Kroupa I J2007l) . we therefore 
de fine Xc.F.xp = 3tm (i.e. we mul tiply the 'gas expulsion time-scale' 
of IBaumgardt & Kroupa I d2007h by a factor 3) so that XcExp corre- 
sponds to a residual gas mass fraction of = 0.05, i.e. the cluster 
is practically devoid of gas. 

We subject all cluster-forming cores to the same external 
tidal field, which is that of an isothermal potential with a circu- 
lar velocity V c = 220km. s^ 1 at galactocentric distances of either 
D ga [ = %kpc or D„ a i = 4kpc. This will allow us to assess how the 
strength of the external tidal field influences modelling outputs. The 
embedded-cluster tidal radius obeys: 



D 



2m ga i(<D 



gal) 



1/3 



(4) 



where m ec i = SFE.m core is the embedded cluster stellar mass and 
ntgai ( < Pool) = Vr Deni /G is the host galaxy mass enclosed within 



D ga i (Bi nnev & Tremaine 111994) . G is the gravitational constant. 

Tidal overflow sets in if the embedded-cluster radius r ec j 
is larger than the tidal radius r,. Substituting m ga i{< D ga /) with 
Vc.Dg a i/G in Eq.[4] it follows that this equates with a volume den- 
sity smaller than: 



Plb 



2kG D 2 , 

gal 



(5) 



At D 



gal 



Skpc, Eq. [5] gives p, in 



6.10~ 24 g.cm- 3 



O.OSMQ.pc , equivalent to an H2 number density 



At D 



gal 



4kpc, these figures are 

3™- ,. _^™-3 



iH 2 Jim = \-1cm 

Plim = 24.10~ 24 g.cm~ : ' = 0.32MQ.pc~3 or riH 2 ji m = 5cm 

Cluster gaseous progenitors are denser than these limits by 
several orders of magnitude. Figure[TJshows mass-radius diagrams 
of molecular cores mapped with different tracers. The top panel 
shows radii and masses of cores mapped with the C 18 7= 1- 
emission line, some of them displaying signs of star formation 
(see Section [3] for a detailed discussion). In contrast, the middle 
and bottom panels present mass-radius diagrams of molecular cores 
selected for their star formation activity, then mapped with higher 
density tracers: CS 7 = 5 — 4 emission line and dust continuum 
emission. We will discuss these observations in greater detail in 
Section[5] For now, suffice is to say that, owing to their systematic 
star formation activity, molecular cores of the middle and bottom 
panels constitute better proxy of cluster gaseous progenitors than 
the C 18 cores of the top panel. Each panel also shows lines of 
constant volume number density («//2> dashed black lines) and of 
constant surface density (£, dotted black lines) fitting the data. The 
mean number density ranges from njj2 = 3.10 cm' - 3 fortheC 18 
cores (top panel), to iihi = 2.10 4 c/n~ 3 (middle panel) and n#2 = 
2. 10 5 cm~ 3 (bottom panel) for the molecular cores selected for their 
star formation activity. Observed molecular cores are thus denser 
than the tidal limit defined by Eq.[5]by about 4 orders of magnitude, 
i.e. they are 'immune' to galactic tides. 

Following gas expulsion, however, gas-loss, infant weight-loss 
and spatial expansion decrease the density of clusters compared 
to that of their parent cores. The key-point our simulations aim 
to address is: in terms of cluster-forming core mass-radius rela- 
tions, what conditions lead to tidal overflow for the expanded clus- 
ters and, therefore, to enhancement of infant weight-loss/mortality 
compared to what would be obtained for isolated clusters (i.e. no 
tidal field). 

We test 6 different mass-radius relations: constant core sur- 
face density, constant core volume density, and constant core ra- 
dius, each with two different normalizations. We parametrize the 
mass-radius relation by its slope S and normalization %: 



1 pc 



X 







(6) 



Table[TJshows adopted % and S values, along with the correspond- 
ing surface densities, volume densities and radii. Models with con- 
stant core surface density (5 = 1/2), constant core volume density 
(5 = 1/3) and constant core radius (5 = 0) are labelled X core , Pcore 
and r core , respectively. For each slope 8, we consider two normal- 
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Table 1. Adopted mass-radius relations r core \pc\ = x{m core [MQ]) s for 
cluster-forming cores, with their corresponding constant surface density 
^■core (5 = 1 /2), volume density p core (5 = 1 /3), and radius r COTe (5 = 0). 



^core 


Pcore 


' core 


5 = 1/2 


5 = 1/3 


5 = 


Compact 0.5g.cm~ 2 


3.10-"g.cm- 3 " 


0.3pc 


x=om 


X = 0.04 


X = 0.30 


Loose 0.05g.cm~ 2 


l(T 20 g.cm- 3 b 


1.5pc 


X = 0.04 


* = 0.11 


Z = 1.50 


" n m = 6-10 4 cm" 3 = p core 


= 41OOM .pc- 3 




b n H 2 = 3. 10 3 cm ~ 3 = p m „ 


= 2OOM .pc~ 3 





izations x, referred to as 'compact' or 'loose' model. The 'loose' 
T-core and p CO re models are fits to the C 18 data with the slope 8 
imposed (dotted and dashed black lines in top panel of Fig. Q}. The 
'compact' Z eore and p CO re relations describe the data of molecular 
cores selected for their star formation activity, for which we adopt 
Ecore = 0.5g.cm~ 2 and nui.core = 6.10 4 cm~ 3 (dotted and dashed 
blue lines with filled circles in middle and bottom panels, respec- 
tively). These densities are at the logarithmic mid-points between 
the data-fits of the middle and bottom panels of Fig.Q] As for the 
constant radius models, we adopt r core = 0.3pc and r core = l.Spc. 
These are shown as the blue ('compact' model, middle panel) and 
black ('loose' model, top panel) solid lines in Fig.Q] In all forth- 
coming fie ures, X core , Pcore and T core models are depicted by dotted, 
dashed, and solid lines, respectively. 

We link the cluster-forming core radius r core and the embed- 
ded cluster half-mass radius r/, with: 



Kr a 



(7) 



We adopt K = 0.5. Molecular cores have power- law density profiles 
p(r) oc r~ p , with th e density index 1 .5 < p < 2. ( IMuller et al.l2002l : 
iBeuther et al J2002h . K = 0.5 corresponds to a truncated isothermal 
sphere (p = 2). Shallower molecular cores (i.e. smaller density in- 
dices p) lead to larger K values since they have a greater fraction of 
their mass in their outer layers. Larger K values in turn lead to larger 
r^/r t tidal field impacts and cores more sensitive to tidal overflow. 
Combining Eqs. [4] [6] and [7] provides the tidal field impact 



n 



0.36 



-2/3 



1M© 
V, 



5-1/3 



220km.s 



2/3 



(8) 



thereby highlighting the influence of both the slope S and normal- 
ization % of the core mass-radius relation. Figure [2] depicts Eq. [8] 
for the various parameters (x, 8) of Table Q] for D ga [ = Skpc (top 
panel) and £> ?fl / = Akpc (bottom panel). For constant core radii 

(8 = 0), more massive clusters sit more deeply within their tidal 

— 1/3 

radii (r/,/r f oc m J ) and are thus more resilient to the external 
tidal field. Conversely, constant surface density cores (5 = 1 /2) are 
conducive to more massive clusters being more prone to tidal over- 
flow ixhl r t x m ld)- I n case °f constant volume density (8 = 1/3), 
the tidal field impact is independent of the embedded-cluster mass 
( r h/ r t °= ra f C /)- Equation[8]can actually be rewritten as a function of 
the core number density njj2 as the sole core parameter. The core 
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Figure 1. Radius vs. mass of dense molecular gas cores. Top panel: molecu- 
lar cores mapped in C ls O emission (filled symbols). Surrounded filled sym- 
bols represent cores with detecte d star formation activit y . Horizontal arrows 
dep ict the lowest-s i ze cor es in lAovama et all 1200 ll) , ISaito et al.1 i 19991) 
and lYonekura et all 120051) , based on their half-power beam width. Middle 
panel : molecular cores map ped in dust continuum emission iMuller et alj 
2002; Faundez et al. 2004(). Bottom panel: molecula r cores mapped in 
CS I Shirley et al. 2003) and dust continuum emission iMuller et al]|2002l ; 
iFontani et al.l l2005t) . Our loose r corc = 1.5 pc, p core = 200MQ.pc~ 3 = 
n H2 = 3.10 3 cm~ 3 and Y, core = 0.05g.cm~ 2 models are shown in the top 
panel as solid, dashed and dotted lines, respectively. The compact mod- 
els r com = 0.3 pc and X. core = 0.5 g.cm~ 2 are shown as the (blue) solid 
and dotted lines with filled circles in the middle panel. The p core = 
4100MQ.pc~ 3 = rifji = 6. 10 4 cm -3 compact model is depicted as the 
dashed (blue) line with filled circles in the bottom panel. 



radius, mass and number density are related by: 

rcore[pc\ = 1.5 — 1 

\n H2 [cm J J 

Inserting Eq.[6]and Eq.|9]in Eq.(8] we obtain 



(9) 



1 
n 



■ 4.2 K SFE 



4/3 -U3fD gal y 2/3 f V, 



l H2 



\ \kpc 



220km. s 



2/3 



(10) 

Higher rjjr, ratios, and thus greater vulnerability to the tidal field, 
result either from a lower core density nf]2, thus larger x ( see loose 
vs. compact models in top panel of Fig. [2}, or from a stronger tidal 
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Figure 2. Relation between the ratio /'/,/(> of the half-mass and tidal radii 
of embedded clusters and the mass m^fire of their progenitor gas core, t 'core* 
Pcore and £ core refer to cores with constant radius, constant volume den- 
sity and constant surface density, respectively. The top and bottom panels 
consider cores exposed to the tidal field of an isothermal halo with a cir- 
cular velocity V c = 220 km. s~ l at galactocentric distances of 8 and 4 kpc, 
respectively. As the loose models (open symbols/black curves) have a den- 
sity lower than that of realistic cluster-forming cores, they are shown for 
illustrative purposes in the top panel only. See text for additional details. 



field, equivalent here to a smaller galactocentric distance (see com- 
pact models in top and bottom p anels of Fig.[2]l. 

Building on Fig. [2] and on iBaumgardt & Kroupa I d2007t) N- 
body model grid of clusters, which provides the fraction F„ oun( i of 
stars bound to a cluster at the end of violent relaxation as a function 
of SFE, TcExp 1 1 Across and r^/rt, we obtain in Fig. [5] the relation 
between F DOUni i and m core . Model parameters are identical to those 
in Fig. [2] and the gas-expulsion time-scale is Xgexd = Across- Note 
the correlation between a low F oollnl j and a high r/Jr, in Fig. [2] 
In Fig. [3] lower normalizations % or larger galactocentric distances 
D qa i result in larger bound fractions through a smaller tidal field 
impact. The bound fraction as a function of mass is constant when 
5 = 1/3, increases when 5 = and decreases when 5 = 1/2. The 
latter illustrates - for the first time - a case where violent relaxation 
preferentially destroys high mass clusters. 

As quoted earlier in this section, the compact mass-radius re- 
lations constitute a better proxy of cluster initial conditions than 
their loose counterparts since they fit data of molecular cores se- 
lected for their star formation activity (see Section [3] for details). 
Therefore, from now on, we focus most of our attention onto the 
compact models. When = 8 kpc, the core mass-radius relation 
and the tidal field impact are essentially two disconnected issues 

Up to Mcore — 

10 5 M o .That is, regardless of the adopted 'compact' 
model (r core , p core or ~L core ), r h /r, < 0.05 when m core < 10 5 M Q . 
In this regime, the tidal field impact is weak and exposed clusters 
respond to the loss of their residual star-forming gas essentially as 



Figure 3. Bound fraction of stars after cluster violent relaxation as a func- 
tion of cluster-forming core mass for the model parameters adopted in Fig. [2] 
and a gas expulsion time-scale set to a crossing-time of the cluster-forming 
core, i.e. toExp = Across- Line- and symbol-codings as in Fig. [2] The thin 
dash-dotted horizontal (black) line in each panel indicates the bound frac- 
tion for the adopted SFE and tGExp / Across in the absence of an external tidal 
field, i.e. F hmmi = 0.35. 



if there were no external tidal field l lBaumgardt & Kroupa II2007T ). 
As a result, the bound fraction is almost constant and indepen- 
dent of the adopted core mass-radius relation (Fbound — 0.3 when 
in-core < 10 Mq, see top panel of Fig. [3}- We remind the reader 
that the constancy of F DOUnc i in this regime also stems from our hy- 
potheses of constant SFE and constant icExpl Across (see Eq.O. At 
masses higher than 10 5 Mq, however, the S fore model on the one 
hand, and the p CO re and r CO re models on the other hand, show very 
different behaviours, with r/,/?> ratios and bound fractions Fj mmi ] 
increasingly different as the core mass increases. A smaller galac- 
tocentric distance (e.g. D„ a i = Akpc instead of D„ a i = Skpc) in- 
creases further the contrast between the p core and £ core models at 
high mass. 

That the p CO re and E eore models show such contrasting be- 
haviours in the high mass regime in a strong tidal field, i.e. close 
to the galactic centre, demonstrates the importance of distinguish- 
ing between these two mass-radius relations. Actually, many spiral 
galaxies show a transition from being predominantly atomic in their 
outer regions to bein g predominantly molecular at their centres 
dWong& Blitz l2002h . One may thus expect that closer to the galac- 
tic centre, the amount of dense molecular gas available to star for- 
mation is larger, thus implying that the cluster-forming core mass 
func tion is sampled up to a higher mas s (size-of-sample effect; see 
also lWeidner. Kroupa & Larsenll2004r) . This in turn would lead to 
the formation of more massive embedded-clusters in stronger tidal- 
field-environments, that is, the regime where the Pcore and H C ore 
models lead to highly different final bound fractions of stars. 

Figure [4] is the mass-radius diagram of molecular cores from 
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Figure 4. Mass-radius diagram for the dense molecular cores of the mid- 
dle and bottom panels of Fig. [T] i.e. those selected for their star formation 
activity. Symbol-coding is as in Fig.[T] The dotted, dashed and solid (blue) 
lines are the adopted 'compact' T. core , p C ore and r core models, respectively. 
The dash-dotted (black) lines are characterised by r/Jr, = 0.15,0.10,0.05 



(from top to bottom) for V c = 220 km.s 
D ga i = 4 kpc (bottom panel). 



and D ga i = 8 kpc (top panel) and 



the middle and bottom panels of Fig. Q] We superimpose onto 
these data the adopted 'compact' mass-radius relations (blue dot- 
ted, dashed and solid blue lines). We also show lines of con- 
stant r/Jr, (black dash-dotted lines) for D ga i = 8 kpc (top panel) 
and D ga i = 4 kpc (bottom panel): r h /r t = 0.15,0.10,0.05 (from 
top to bottom). Note that iso-r/,/r> lines are vertically shifted by 
Alogr core = —0.2 in the bottom panel compared to the top one 
stnce r h /r, « Z^f (see Eq.®. Note also that i S0 -r h /r t lines are 
lines of constant volume density, as shown by Eq.fTOl 

For the galactocentric distances considered, the vast major- 
ity of the observed molecular cores are tidal-field resilient, i.e. 
r h/ r t < 0.05. The observational data, however, occupy a limited 
mass range, with only a few cores more massive than 3.1O 4 M0. 
They therefore fail to probe the high-mass regime where we predict 
Pcore and ^eore models to respond very differently to gas expulsion 
through the tidal field. 

Figure [4] allows us to understand Figs. [2] and [3] from another 
perspective. While the compact p core model has r/Jr, < 0.05 irre- 
spective of core mass, a L core model is characterised by a decreas- 
ing mean volume density p core with increasing core mass. This 
equates with a greater ri,/r, tidal field impact for more massive 
cores (Eq. [10] and Fig. [4]l and, thus, a lower final bound fraction 
of stars Ftound (Fig.[3j. Conversely, cores with mass-independent 
radii (r core model) increase their volume density with their mass, 
rendering less massive cores more prone to tidal overflow through 
a larger ri t /r t ratio. 

One may argue that the reason why cluster infant weight- 
loss/mortality for the E core model is so prominently mass- 
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Figure 5. Bound fraction of stars at the end of violent relaxation as a func- 
tion of the core mass. Top and middle panels: Same as Fig. [3] for slower 
gas expulsion (Tq Exp / T cross = 3). Bottom panel: Same as top panel but for 
SFE = 0.25. 

dependent in Fig.[3]partly stems from the high adopted upper limit 
on the core mass range, i.e. log(m cori? [MQ]) = 7.5. One should keep 
in mind, however, that such a large mass of dense molecular gas is 
needed to form a cluster of a few million solar masses, that is, with 
a mass comparable to that of the most massive old globular clus- 
ters and star clusters formed in galaxy mergers. For the p core model 
in Fig. [3] SFE ~ 0.3 and Fj, ounc i ~ 0.3 lead to a cluster mass m c / 
at the end of violent relaxation an order of magnitude lower than 
its progenitor core mass (see Eq. [TJ. As for the 'compact' ~L core 
model at D ga i = 4kpc, it prevents the formation of massive clus- 
ters since cores more massive than ~ 10 6 Mq fail at forming bound 
clusters (Fi, ounc [=0), and 3. 10 5 Mq cores give rise to bound clusters 
~ 10 4 Mq in mass only (SFE ~ 0.3 and Fi, olmc [ ~ 0.1). 



2.2 Exploring a wider parameter space 

Figure [5] explores more widely the parameter space (SFE, 
tGExp/ Across)- Its top and middle panels are the counterparts of 
Fig.[3] with longer gas expulsion time-scales: tQExp — Across- For 
slower gas expulsion, clusters retain a higher fraction Ft, olmc i of their 
stars because they are better able to adjust to the new gas-depleted 
potential they sit in. Besides, slower gas expulsion is conducive 
to smaller spatial expan sion of the exposed cluster (see fig. 3 in 
iGever & Burkert 11200 lh . thus to a higher mean density at the end 
of violent relaxation and greater resilience to the external tidal field. 
This is another channel through which the bound fraction of stars is 
increased compared to quicker gas expulsion. Compared to Fig. [3] 
the bound fractions Fi, mmc i in top and middle panels of Fig.[5]are in- 
creased by factors 2.5 — 3. At a galactocentric distance £) ?a / = Skpc 
(top panel), this strongly dampens any dependence of the bound 
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Figure 6. Top and middle panels: Final bound fraction of stars as a func- 
tion of the core mass for the compact (top) and loose models (middle) in 
case of long gas expulsion (TaExp — 30t [TOJ! ). Symbol-coding as previously: 



models, triangles; p core models, circles; 



models, squares. Bot- 



tom panel: Mass-radius relations for the loose and compact L core and r core 
models. The thick (black) dash-dotted line is the minimum number density 
molecular cores must have for their clusters not to be affected by tides in 
case of adiabatic gas expulsion with SFE=0.25 and D ga / = 8 kpc. The ver- 
tical arrows actually indicate that the crossing of that density limit in the 
mass-radius diagram corresponds to a sharp increase of the bound fraction 
of stars. 



fraction on the core mass for the £ core model (0.5 ^ Fbound ^ 0.75). 
But it also results in cluster infant weight-loss hardly compati- 
ble with ob servations since inf ant weigh t-loss is reported to range 
from 70 % dBastian et ai]|2005h to 90 % jLada & Lada1l2003l) . that 
is, 0.1 ^ Fbound ^ 0.3. The bottom panel of Fig. [5] illustrates the 
Fbound vs - m core relation for a lower SFE, namely, SFE = 0.25 
and the same gas-expulsion time-scale. In that case, cluster sur- 
vival {Fbound > 0) requires TGExp > Across- We will further discuss 
the consequences of these plots for the cluster mass function in 
Section |4~T1 Note that the models of the middle and bottom pan- 
els behave almost similarly, that is, the combination of the weaker 
tidal field and lower SFE in the bottom panel compared to the mid- 
dle one leads to a model degeneracy. 



2.3 Adiabatic gas expulsion 

Adiabatic gas expulsion (tcExp » Across) allows an analytic anal- 
ysis which illuminates these results further. The top and middle 
panels of Fig. [6] depict the evolution of Fbound with tne core mass 
for the compact and loose models, respectively, for D ?n / = 8 kpc, 



SFE = 0.25 and TGEXp/ Across = 3 (i.e. the longest gas expul - 
sion time-scale in the model grid of iBaumgardt & Kroupa 112007b . 
Line- and symbol-codings are as in Figs. [5] and [5] For isolated 
clusters (r/Jr, < 0.05), so long a gas expulsion time-scale when 
SFE = 0.25 is conducive to Fb oum { — 1, namely, no cluster infant 
weight-loss. In fact, in a tidal- field-free enviro nment, adiabatic gas 
expulsion implies Fb oun d — 1 < lMathieu1ll983l) . If there is a strong 
enough tidal field, however, stars driven beyond the cluster tidal 
radius by its spatial expansion will get unbound and Fb ou „d < 1. 
In the framework of the adiabatic approximation, we can estimate 
what minimum number density njj2 the cluster progenitor core 
must have to prevent tidal overflow. 

In the case of adiabatic gas expulsion, the radius multiplied by 
the mass is an adiabatic invariant: the cluster expands by a factor 
SFE -1 after gas expuls ion while its mass gets lowe r than the core 
mass by a factor SFE faills II 198tj ; iMathieu |[l983l) . The gas-free 
cluster density thus follows: 



3 m c [ 



SFE.m c 



4iz (SFE' 



•'core ) 



-SFE\p c , 



(11) 



Therefore, for SFE=0.25 and adiabatic gas expulsion, an expanded 
cluster at D ga i = Skpc experiences tidally-driven mass-loss (p c / ^ 
Plimit ' Eqf5J if its parent core has p core < SFE 4 x p; iffl , or n Hl < 
300cm~ 3 . In contrast, clusters formed out of cores with «// 2 > 
300 cm~ 3 are not or little affected by tides. 

This density limit is shown as the thick dash-dotted (black) 
line in the r core vs. m core diagram of Fig. [6] (bottom panel). In this 
diagram, the intersections of the nn 2 = 300 cm -3 density limit with 
lines of constant T. core give the core mass above which constant 
surface density cores give rise to clusters significantly affected by 
tides. Similarly, the intersections with lines of constant r core gives 
the core mass below which clusters formed out of constant radius 
cores experience tidal overflow. These intersections are highlighted 
by upside-down open triangles in the bottom panel of Fig. [6] 

Let us consider the loose model Y. core = 0.05 g.cmT 1 . Its 



intersection in the m„ 

,-3 



- r core diagram with the density limit 
riH 2 = 300cm - - 3 renders log(m eore ) ~ 4.7. This matches the 
core mass regime over which Fb oun d drops significantly for that 
particular core mass-radius relation (open triangles in the middle 
panel of Fig. |6), as indicated by the vertical dotted double-head 
arrow. log(ra co , e ) < 4.2 corresponds to n// 2 > 300cm -3 for which 
we expect little or no infant weight-loss. The middle panel of 
Fig. [6] indeed shows that 0.75 < Fb oun d < LOO over that mass 
regime. In contrast, log(m cnre ) > 5.2 leads to nn 2 < 300cm~ 3 , 
for which we expect expanded clusters to be severely affected by 
tides. The middle panel of Fig. [6] indeed predicts Fbound — for 
log(m core ) > 5.2. A stronger tidal field (i.e. closer to the galactic 
centre) would lower the mass-limit at which Fbound decreases. 
Similarly, the intersection between the constant radius loose 
model r core = 1.5 pc and >ih 2 = 300 cm yields log(/?/ core ) — 2.5, 
where Fbound ls sharply increasing from to 0.75, as indicated 
by the solid double-head arrow. As for the compact models, 
which better describe star cluster initial conditions, those are more 
resilient to tidally-driven mass-loss. The m core — r core diagram 
confirms that Fb oun d must strongly decrease at log(ra corf ) ~ 7.5 for 
^core = 0.5 g.ctrT 2 (see dotted upward arrows in top and bottom 
panels of Fig.|6]l. 

We underline that the above-described effects take place even 
if all cluster-forming cores are located within a galaxy region over 
which the external tidal field does not vary markedly. In these mod- 



8 G. Parmentier & P. Kroupa 



els, the Ft, ounc i(m core ) variations are solely driven by the core mass- 
radius relation. 

Although unbound stars located beyond the cluster tidal radius 
linger around the cluster and may result in observed clusters with 
estimated radii larger than their tidal limit, observations of young 
clusters in the Small Magellanic Cloud show that this effect fades 
away by an age of about 30 Myr dGlatt et al.l2010h . 



3 WHAT INFORMATION CAN WE EXTRACT FROM 
THE OBSERVATIONS OF DENSE MOLECULAR GAS 
REGIONS? 

In Fig. [TJ we compile masses and radii of dense gas regions from 
the literature. The top panel shows results of C 18 mapping of 
dense gas regions in neighbouring giant molecular clouds (GMC). 
The middle and bottom panels are mass-radius diagrams of molec- 
ular cores selected for their star formation activity. 

Thre e Galactic regions are encompassed by the top panel: 
Orion B l lAovama et aiTl200lh . the GMC to ward HII regions S35 
and 37 (referred to S 35/37 in what follows) dSaito et alj|l999l) and 
the Tj Carinae GMC dYonekura et alj|2~005l) . at assumed distances 
of 400 pc, 1.8 kpc and 2.5 kpc, respectively. All observations were 
performed at the same resolution with the NANTEN telescope. 
The lower bound on resolved-core radii for each data-set as im- 
posed by the NANTEN half-power beam width (2'. 7) is shown by 
dotted horizontal arrows. That C 18 cores in the r\ Carinae GMC 
are larger than in Orion B is thus a purely resolution-driven effect. 
C 18 emission traces molecular gas with number densities nu 2 = a 
few 10 3 - 10 4 cm , and the volume density range covered by these 
observations is thus an imprint of the C I8 tracer. 

The 

Score an d Pcore loose models of Section [2] are fits to 
these C 18 data. As we cautioned in Section[2] C 18 cores do not 
systematically host signs of star formation activity and, therefore, 
loose models may not trace actual cluster formation conditions. In 
Fig.[T] C ls O cores associated to one or more IRAS sources iden- 
tified a s a Young Ste l lar Ob je ct (YSO) candidate by ISaito et all 
dl999h , lAovama et alJdioOlh or lYonekura et al.l ( 120051) are circled. 
Open squares depict 77 Carinae cores showing other signs of active 
star formation (e.g. YSO candidate from the MSX point-source cat- 
alog or bipolar outflows). 

c 18 o cores were also observed by I Aovama et al 

]|IIl|) 

and 

lYonekura et all fcOOSl) in H 13 CO+ 7=1-0 emission. H 13 CO+ 
emission traces molecular gas at njj 2 — 10 5 c/n~ 3 , i.e. an order of 
magnitude higher than C 18 0. They note a tight correlation between 
the presence of high-density H 13 CO + clumps and star formation 
activity. This suggests that star formation requires number densities 
of order nn 2 — 10 5 cm~ 3 (= p CO re — 7000 Mq ./jc~ 3 ) and hence 
that only the densest, presumably most inner, regions of C ls O 
cores form stars. Loose Z C ore and p core models have therefore 
too large a normalization x to emulate realistic cluster formation 
conditions. That is why we insisted in Section[2]that those models 
should be considered only for illustrative purposes, e.g. to show 
how model outputs respond to variations of the normalization 
of the core mass-radius relation. A related point worth being 
quoted here is th at SFEs measured o ver the whole volume of 
C 18 cores (e.g. Illiguchi et alj I2009L their table 3) are global 
SFEs and, as such, are not indicative of how an embedded cluster 
dynamically responds to gas expulsion. This is the local SFE, 
namely, the SFE estimated over the volume of gas forming the 
cluster, which matters when modelling cluster violent relaxation. 
Global SFEs averaged over whole C 18 cores constitute lower 



limits to their local counterparts. A low global SFE (say, 10 %) 
may be misleading in prompting us to conclude that a cluster will 
not survive its violent relaxation, even though the local SFE may 
be high enough for the cluster to retain a fraction of its stars. We 
will com e back to this point in a forthcoming paper jParmentier. I 
lin prepl) . 

To better constrain cluster formation conditions, we gather in 
the middle and bottom panels of Fig. Q] masses and radii of dense 
molecular cores selected for their star formation activity (either 
IRAS sources or water masers). Mapping of star-for ming cores in 
the CS J= 5 — 4 emission line has been performed bv lShirlev et al.l 
( I2003L columns 3 and 5 of their table 5; filled squares in the bottom 
panel of Fig. [T}. Mapping of st ar-forming cores in du st-continuum 
emission has been p erformed by lFaundez et ai]d2004l their table 1), 
iFontani et al. 

1( 120(4 

radii and masses from their tables A.5 and A.6, 
respectively) and iMuller et alj ( 120021 columns 2 and 3 of their ta- 
ble 4). They are depicted as the (blue) x -symbols and asterisks and 
(black) open circles in the middle and bottom panels of Fig. Q] The 
radius of cores is defined as that of the contour at half-maximum of 
the CS or dust-continuum emission. The core mass is the mass en- 
closed within that radius. These cores have (volume and surface) 
densities significantly higher than those inferred from the C 18 
data, as indicated by the lines of constant E C ore an d constant nu2 
in middle and bottom panels. iMtiller et alj d2002l) provide an al- 
ternative definition of the masses and radii of their surveyed cores 
(columns 4 and 5 of their table 4), which we show as the (black) 
filled circles in the middle panel of Fig. Q] This r core — m core se- 
quence neatly defines a line of constant volume density, with a 
mean number density n# 2 ~ 2.10 4 cm~ 3 . 

It is important to realize that this result stems from how core 
masses and radii are estimated, however. From th e radial den- 
sity pr ofiles of the star-forming regions they observe, iMiiller et alj 
j2002l) obtain the radius where = 10 4 cm 3 . Their core mass is 
the gas mass enclosed within that radius. Defining core masses and 
radii that way necessarily results in a m core — r core sequence of con- 
stant volume density. As such, this core mass-radius relation may 
constitute a measurement-imprint rather than a genuine imprint of 
the cluster-formation physics. 

These various examples show that to infer the mass-radius re- 
lation of cluster-forming cores observationally is not a straightfor- 
ward task. Results heavily depend on the tracer and/or the method 
used to map them. In Sections|2]and|4]realistic cluster initial condi- 
tions are described by mass-radius relations representative of the 
dense molecular cores selected for their star formation activity. 
We refer to them as the 'compact' £ C ore, r core an d Pcore mod- 
els. They are shown as the (blue) dotted, solid and dashed lines 
with filled-circles in the middle and bottom panels of Fig.Q] Their 
constant surface density, radius and volume number density are 
Score = 0.5 g.cm~ 2 , r core = 0.3 pc and n//2,core = 6 x 10 4 cm~ 3 , re- 
spectively. These volume and surface densities are at the logarith- 
mic midpoints of the fits to the data in Fig. Q] middle and bottom 
panels (black dashed and dotted lines). 

We emphasize that our core mass-radius relations (Eq. [6] 
and Table QJ relate the total mass of cores to their outer ra- 
dius. In that sense, our relations are not di rectly comparable t o 
the dust-continuum and CS emission d ata of IShirlev et alj d2003l) . 
iFaundez et alj d2004l) and lFontani et alj d2005l) . who define the core 
as the region enclosed within the FWHM contour. Assuming an 
isothermal sphere density profile for the cores, the volume and sur- 
face densities within the half-mass radius are 4 and 2 times, re- 
spectively, higher than the volume and surface densities averaged 
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over the whole core (since r/, = 0.5r core for an isothermal sphere). 
That is, the mean volume and surface densities within the half- 
mass radius of our compact E core and p CO re models are fully com- 
parable to the mean densities within the FWHM contour of the 
data of Fig. [Tj bottom panel. Note also that the adopted number 
density «H2.eore = 6 x 10 cm -3 is not significantly different from 
the number density characterizing H I3 CO + -traced molecular gas 
( n H2,core — 10 5 cm -3 ) which is closely associated to star formatio n 
activity in C ls O cores jAovama et al.ll200lllYonekura et alj2005l) . 

Physically, a mean density of n# 2 ~ 10 4 -10 5 cm~ 3 (~ 700- 
7000Mq./)c~ 3 ) for cluster formation may result from the asso- 
ciated efficient decay of turbulence, leading such dense gas cores 
to un dergo gravitational collapse and form star clusters dKlessen I 
120031) . We also note that n Hl ~ 10 5 cm -3 leads to z cmss ~ 0. \5Myr, 
implying that all cluster stars due to become unbound owing to gas 
expulsion have crossed t he tidal radius bou ndary by an age of at 
most 15 Myr (see fig. 4 in lParmentier 112003 . see also Section POl l. 



4 FROM CLUSTER EARLY EVOLUTION TO GALAXY 
STAR FORMATION HISTORIES: CONSEQUENCES 

Section [2] shows that the combination of the cluster-forming core 
mass-radius relation with an external tidal field can contribute to 
determining how much infant weight-loss clusters experience and 
whether infant weight-loss is mass-independent or not. In this sec- 
tion, we survey a few topics which are directly influenced by the 
core mass-radius relation, either in relation to the tidal field impact 
(Sections |4. 1 | and [4~2t . or independently of it (Sections |4.3| and |4~4l i. 

4.1 The shape of the young cluster mass function 

Most observational evidence gathered so far shows that the shape 
of the post-violent relaxation cluster mass function mirrors that 
of the embedded c l uster mass function jKennicutt et al.l 19891; 
McKee & Williams I [l997l: lLada & Ladal 120031: IZhang & Fall! 



19991 : lOev et all |2004 iDowell et all 120081 , but see Anders et 



al. (2007) for the case of a bell-shaped young cluster luminosity 
function). That is, cluster infant weight-loss appears to be mass- 
independent. As Figs.[3]and[5]show, under the assumption of con- 
stant, hence mass-independent, SFE and gas expulsion time-scale 
iGExpl Across, cluster-forming cores characterised by a constant vol- 
ume density (5=1 /3) constitute the most robust way of achieving 
mass-independent cluster infant weight-loss. In contrast, in the case 
of constant surface density, clusters formed out of massive cores 
may be preferentially destroyed (see Fig. [3j, and the shapes of the 
post-violent-relaxation cluster mass func tion and core mas s func- 
tion may differ substantially (see fig. 4 in lParmentier Il20ich . 

Our result is at odds with that derived by 
iFall. Krumholz & Matznerl j201Ch . following which mass- 
independent infant weight-loss requests near-constant surface 
density cores (5 ~ 0.5). In the case of constant volume density 
cores, they find that the SFE needed to clean the cluster of its 
residual star-forming gas is an increasing function of the core 
mass. This is conducive to less-massive clusters experiencing 
greater infant weight-loss and, thus, to a cluster mass function 
shallower than the embedded-cluster mass function. Our result and 
theirs stem from two utterly different approaches, however. Our 
model rests on how the tidal field impact r^/r t varies with the core 
mass m core , under the assumptions of constant SFE and constant 
gas expulsion time-scale tQExpl Across- Their model rests on the 
amount of stellar feedback required to clear an embedded cluster 



of its residual gas, neglecting the tidal field impact and assuming 
constant ZoExp/^cross- 

To illustrate that both approaches are not irreconcilable, 
let us first consider the case of energy-driven feedback of 
IFall. Krumholz & Matzner 1 d2010l) . The rate E at which massive 
stars deposit energy in the cluster-forming core gas is proportional 
to the core stellar mass, that is, E = kE-SFE.m core with Ice a pro- 
portionality coefficient. The energy input accumulated over the 
gas expulsion time-scale IpExn ls t hus E tot = kg.SFE.m core .%QE xp . 
IFall. Krumholz & Matzner derive the core SFE by equat- 

ing the total energy input E, ot to the critical value needed to ex- 
pel the intra-cluster gas, that is, the gas binding energy E cr j t ~ 
G(l — SFE) m^ ore /r core . Introducing the core crossing-time t cmss , 
it thus follows: 



Eiot 

Fcrit 



1 cz7i- ^GExp 
■ k E ■ SFE . mcore ■ ■ T« 



G(i-SFE)m, 



Across 
2 

core 



r core 



(12) 



Since T crojs = k r rclf e /(Gm core ) 1 / 2 , with k r a unit-dependent pro- 
portionality constant: 



(k E .k T .G-^ 2 ). 



SFE TqExp 



I -SFE TV, 
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: m core ■ f core ■ 



(13) 



Finally, introducing the core mass-radius relation (Eq.[6}: 



( X 5/2 .k E .k t .G 



-3/2n 



SFE Zgex 



1 SFE Tcross 

_ m (3-5«)/2 
— " l core 



(14) 



Neglecting the coefficient 1 — SFE which matters little as long as 
SFE < 0.5, we see that, depending on the slope 5 of the core mass- 
radius relation, the product SF 'E .ZgExp I Across increases with the 

core mass as mJcte (5 = 1 /2), rriJore (5 = 1 /3), and rr?J ^e (<5 = 0). 
The greater dependence of SFE .Xqe xp / Across on the core mass as 
5 decreases stems from the depth of the core potential well being 
itself a steeper function of the core mass for shallower core mass- 
radius relations. 

Equation[l4]can now be applied to two limiting cases: either a 
mass-independent SFE, or a mass-independent gas expulsion time- 
scale in units of a core crossing-time. A c onstant SFE is the ap- 
proach adopted bv lBaumgardt et al ] J20081) . which we come back 
to below. 

If, on the other hand, Tqexp I Across ls core-mass-independent, 
Eq.[l4]leads to 



SFE oc m c lre 5S ^ 2 



(15) 



and w e have recovered eq. (la) of IFall. Krumholz & Matznerl 
Their result suggests that mass-independent cluster infant 
weight-loss demands 5 ~ 0.6. That is, compared to cores with con- 
stant volume density and constant radius, constant surface density 
cores introduce the smallest mass-dependence for cluster-infant 
weight-loss. Yet, Eq. [2] shows that the bound fraction F DOUnc i of 
stars at the end of violent relaxation does not depend on SFE and 
tGExp I Across only. It also depends on the tidal field impact which, 
as we demonstrate in Section[2] can introduce a strong core-mass- 
dependence for 5 ~ 0.5. A mass-independent tidal field impact 
ri,/r, requests constant volume density cores, but those are charac- 
terised by a mass-dependent SFE to expel residual star-forming gas: 
SFE oc mJore (Eq.ll5t. In that case and for constant gas expulsion 



10 G. Parmentier & P. Kroupa 



3.5 4 



log ( m core [MJ) [Eq. 18] 
4.5 5 5.5 6 6.5 7 



7.5 



0.75 



0.5 



0.25 



r h' r i ■■ 



=0.03 



SFE=0.25 
SFE=0.30 
SFE=0.33 
SFE=0.35 
SFE=0.40 



-2 



-1.5 




lo S ( ^GExp^cross ) 



Figure 7. Bound fraction of stars at the end of violent relaxation as a func- 
tion of the gas expulsion time-scale TgExp I Across expressed in units of a core 
crossing-time for the SFEs quoted in the key and a weak tidal field impact 
rf,/r t = 0.03. Based on the W-body model grid of iBaumgardt & Kroupa I 
120071) . The gas expuls ion time-scale Tqexd is de fined as three times the 
e-folding time %m of Baumgardt & Kroupa I ^0071) (see Eq.fJ) 



time-scale tGExp I Across (as assumed in Eg. 115b. Ff, rmf j is a sharply 
increa sing function of the SFE. Figure 1 in IBaumgardt & Kroupa I 
J2007t) shows that Fb„„„d = when SFE < SFE,^ and < Fbound ^ 
1 when SFE, i, < SFE ^ 1, with SFE,/, a threshold value depen- 
dent on icExpl tcwss and r^/r,. For instance, explosive gas expul- 
sion (tcExp « Across) and no tidal field impact (r/Jr, < 0.01) ren- 
ders SFE tn ~ 0.33. The transition from Ff, ounc i = to Ff, ounc i > as 
the SFE increases beyond the threshold SFE,/, is conducive to the 
formation of features in the cluster mass function (flattening and 
turnover), in conflict with most observations of young star clusters 
in the present-day Universe. 

As a brief summary before heading further: under the 
assumption Xqexp 1 1 Across = constant, 8 = 0.6 results in a mass- 
independent SFE (Eq,ll5t. but a mass-dependent tidal field impact 
ri,/r, (Eq.[8}. Conversely, 5 = 1/3 results in a mass-independent 
tidal field impact o»/r f , but a mass-dependent SFE. These 
modelling results are to be contrasted with the observations of 
power-law mass functions for young star clusters which demand 
that all 3 parameters - SFE, TcExpl Across and r^/r, - weakly 
depend on the core mass m core to ensure that the bound fraction 
Fbound does not depend significantly on m core (Eq.[2]l. 



The approach adopted by IBaumgardt eTaD d2008l) can help 
us solve this intriguing conundrum . Instead of assuming a mass- 
independent TQExp/fcross m Eq. 1141 IBaumgardt et alj 12008) adopt 
a mass-independent SFE. Equation [14] thus becomes: 



*GEx 



„(3-55)/2 



(16) 



While in Eq. Q3] the larger energy-input required to clear the gas 
out of more massive cores arises from a higher SFE, in Eq.QJo] it is 
obtained by integrating the energy input over a longer gas expulsion 
time-scale tcExpl Across- Results obtained for TcExp I Across based on 
Eq.[T6]are at first glance similar to those obtained for the SFE based 
on Eq. [15] A mass-independent gas expulsion time-scale requires 
8 ~ 0.6, which leads to mass-dependent tidal field impact. Con- 
stant volume density cores (8 ~ 1/3), needed to reproduce a mass- 
independent tidal field impact, are conducive to mass-dependent 
tGExp/tcross- There is a major difference with the SFE-varying ap- 



proach of Eq.[T5] however. Whether the mass- varying TGExp 1 1 Across 
of Eq. QJo] induces a mass-dependent Fb OU nd depends very much 
on the rang e of Tr.Exn I Across involved. A ctually, the N-ho&y sim- 
ulations of IBaumgardt & Kroupa ] J2007I) show that for XcExp ^ 
0.5r crojJ , the bound fraction of stars stays about constant. That is, 
it is doable to get F hmmd {SFE, %GExpl Across > r h/ r t) = constant even 
though TcExp I Across is an increasing function of the core mass, as 
long as T GExp < 0.5 

Across • 

To assess this issue in more detail, let us derive the normal- 
izing factor in Eq. ll6| Building on mode ls of deposition of stellar 
feedback energy. IBaumgardt et alj ( T2008I) derive the gas expulsion 
time-scale Zqexp ( m units of Myr) as a function of the core half- 
mass radius, core mass and SFE (their eq. 14 which we reproduce 
below for the sake of clarity): 



*GEx 



1.1 x 10 



_o 1 — SFE m cc 



Myr . 



(17) 



SFE M 

EquationQTJcan be combined with the core crossing-time and 
with a core mass-radius relation to infer Ig exb 1 1 Across as a function 
of the sole core mass m core . Using eq. 6 in IBaumgardt & Kroupa I 
d2007h for the core crossing-time and the core mass-radius relation 
of our compact p CO re model (Eq.|6]with 5 = 1/3 and % = 0.04), we 
obtain: 
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(18) 



This equation is valid for constant volume density cores with the 
normalization % = 0.04. 

Figure [7] presents the bound fraction Fi, ounc i in dependence 
of the gas expulsion time-sca le Hqexo I Across (bottom x-axi s, based 
on the /V-body model grid o f IBaumgardt & Kroupa ||2007]) and of 
the core mass m core (top x-axis, based on Eq. 1 1 8b - The adopted 
tidal field impact is weak, namely, r^/r, ~ 0.03, as we find for the 
compact p core model in Fig. [2] One can see that SFE of ~ 0.35- 
0.40 leads to a constant bound fraction up to m core ~ 3 x 1O 5 M0, 
and to an increase by a factor of ~ 4 over the high mass range 
iticore 3. 10 7 Mq . An increase of Fbound by a factor 
of 4 may be hard to detect in the cluster mass function of realis- 
tic cluster systems. Those indeed suffer from the intrinsic scatter 
of cluster individual properties and from measurement uncertain- 
ties. In any case, variations of the bound fraction Fbound by a factor 
of a few constitute a much weaker mass-dependence than a transi- 
tion from a zero- to a non-zero Fb oun d, as yielded by the tidal field 
impact upon constant surface density cores (e.g. bottom panel of 

Fig. ED. 

If our argument that cluster-forming cores have a constant vol- 
ume density, thus mass-independent tidal-field impact, and con- 
stant SFE, is correct, then Fig. [7] suggests that the SFE at the 
onset of gas expulsion cannot be lower than ~ 0.35. Actually, if 
SFE = 0.33, the formation of a bound cluster, able to survive its vi- 
olent relaxation, requires Xqexp > f eTO j s /3, thusow > 3.1O 5 M . 
Lower mass cores, characterised by quicker gas expulsion, fail 
to form a bound cluster, i.e. Fb OU nd = 0. This equates with the 
mass-function of fooMraii-cluster-forming cores to be truncated at 
~ 3.10 5 Mq, which is co nducive to the formation of a turnover in 
the cluster mass functio n jParmentier & Gilmore J2007l) . Similarly, 
IBaumgardt et alj d2008l) and lParmentier et alj d2008l) infer that the 
mass-varying gas expulsion time-scale of constant radius cores can 
result in cluster mass functions substantially different from the em- 
bedded cluster mass function, provided that the SFE is low, namely, 
SFE < 0.35-0.40. 

In a second approach, based on a momentum-driven feedback 
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model, iFall. Krumholz & Matzner their eq. lb) derive the 

SFE such that the momentum imparted to the gas reaches the crit- 
ical value p cl -i, = m core x V e , where V e is the core escape velocity. 
That is, the velocity V s f, e u of the shell collecting the residual star- 
forming gas is the escape velocity V e at the time it reaches the core 
edge. This implies that the gas is driven out on a time-scale of order 
a crossing-time, i.e. Igexd I Acr oss — 1 regardless of the cor e mass . 
Note that equations (la-b) of IFall. Krumholz & Matzner I J2010h 
being written as proportionalities, they are valid when IcExp x 
Across, and not just when %GExp — Across- The energy-driven feed- 
back model, however, shows that gas can be expelled on time-scales 
shorter or longer than a crossing-time (see Fig. |7). One may be 
puzzled by a gas expulsion time-scale Tqexp longer than a core 
crossing-time T cross since it implies a gas velocity slower than the 
core escape velocity V e . Yet, this is not necessarily surprising. The 
escape velocity is defined for objects subjected solely to gravity, 
while the shell of gas, in addition to being subjected to its own self- 
gravity, keeps being powered and pushed outwards by continuing 
stellar winds, ionized region over-pressure, etc. In other words, a 
shell velocity slower than the core escape velocity may not neces- 
sarily preclude the gas from being unbound from its parent core. 



4.2 The star formation history of galaxies inferred from the 
age distribution of their surviving star clusters 

Within the paradigm that most stars form in gas-embedded clus- 
ters, cluster infant weight-loss and infant mortal ity appear as sig- 
nificant drivers of galaxy field star populations jKroupa & Boilv I 
|2002| ; IParmentier & de Griis1l2008l) . Consequently, reconstructing 
the star formation history of galaxies from their cluster age distri- 
bution requires a firm grasp on the time-evolution of the integrated 
bound fraction F BOUru t, that is, the bound fraction F oolm i of stars in- 
tegrated over the core mass function at a given age (see fig. 5 in 
IParmentierl I2009T) . As we now illustrate, this also requires a fair 
knowledge of the core mass-radius relation. 

Building on the compact models of Fig. [3] (i.e. D ga i = %kpc 
and D ga i = Akpc), Table [2] illustrates how F DOlm d varies with 
the slope 8 of the core mass-radius relation, and with the up- 
per limit m up and slope — /3 of a power-law core mass function 
dN oc m~Pdm. We consider two upper limits to the core mass 
range to illustrate the sensitivity of /*/,„,„„/ to this parameter. In the 
left part of Table d m U p = 3 x 10 4 Mq is the upper limit to the 
mass range of the data shown in Fig. Q] and from which we obtain 
the normalizations of our core mass-radius relations. In the right 
part, we adopt m up = 10 7 Mq. While most reported young clus- 
ter mass functions have /3 ~ 2 (see references in Section |4T) , the 
mass function of molecul ar cores is shallower with /3 ~ 1 .7 jLadal 
ll992l ; lKramer et alJI 19981). Given that uncer tainty, Table |2] consid- 
ers both cases (see also lParmentier. lin preph . When j3 < 2, higher- 
mass cores contribute a greater fraction of the total gas mass than 
their low-mass counterparts. For instance, jS = 1.7 leads to cores 
more massive than 10 6 Mq to contribute ~ 50% of the total gas 
mass for a core mass range 10 2 - 10 7 Mq . Under our assumption of 
a constant SFE, this equates with the most massive cores contain- 
ing most of the stellar mass. Therefore, if jS = 1.7, the total stellar 
mass fraction which stays in bound gas-free clusters is sensitively 
determined by the fate of the most massive embedded-clusters. 

For p core models, F oounl j is independent of the core mass 
function parameters /3 and m up since clusters experience mass- 
independent infant weigth-loss. As quoted earlier, the compact fcore 
and p core models do not lead to significantly different F DOlml j(m core ) 



relations, which are also similar for both D c;a i = %kpc and D ga [ = 
4kpc (Fig.[3j. This is a direct consequence of having 7^/r< < 0.05 
(no or weak tidal field impact) over the whole core mass range 10 2 - 
10 7 Mq for the r core and p CO re models when D ga \ > Akpc (Fig. |4j. 
As a result, the integrated bound fraction is Fi, oun d ~ 0.3 for r core 
and Pcore models at D ga [ = 4-8 kpc, irrespective of S and m up . 

In contrast, to raise m up from 1O 4,5 M0 to 10 7 Mq reduces 
Fbound f° r ^eore models since in that case higher-mass embedded 
clusters are more efficiently destroyed than their low-mass coun- 
terparts. Besides, this sensitivity of the integrated bound fraction 
Fbound t0 me core mass upper limit is strengthened for shallow 
power-law core mass functions (j3 = 1.7) and/or in case of stronger 
tidal field (e.g. _D ?a / = Akpc instead of D ?a / = Skpc). 

Independently of the tidal field impact, we note that the 
cluster-forming core mass-radius relation is also relevant to the two 
topics discussed in the next sections. 



4.3 The evolutionary rate of star clusters through violent 
relaxation 

The mass-radius relation of cluster-forming cores also determines 
their mean density and crossing-time T cross , thus how fast clusters 
experience infant weight-loss following gas-expulsion: the shorter 
the core crossing-time, the faster cluster evolution through violent 
relaxation. If cluster-forming cores have a constant surface den- 
sity, more massive clusters evolve more slowly than their low-mass 
counterparts owing to their lower volume density, hence longer 
crossing-time. In contrast, if cluster-forming cores have a constant 
radius, more massive clusters have a higher density, thus shorter 
crossing-time and the duratio n of their violent r elaxation (in units 
of Myr) is shorter (see fig. 1 in lParmentier I d2O10h for an application 
to the time-evolution of the cluster mass function). Constant vol- 
ume density cores are all characterised by the same crossing-time 
and, thus, evolve at the same rate through violent relaxation regard- 
less of their mass. If the mean number density of cluster-forming 
cores is nn 2 c± 10 4 cm , Tcross — 0.45Myr and all cluster stars due 
to become unbound owing to gas expulsion have crossed the tidal 
radius boundary b y a cluster age of at most 100T crojs (see fig. 4 in 
IParmentier 1 120091) or 45 Myr. If the core mean number density is 
10 times higher, ~ 10 5 cm -3 , then Tcross — 0.15Afyr, implying 
that violent relaxation is over by an age of ~ 15 Myr. 



4.4 The mass-metallicity relation of old globular clusters 
predicted by self-enrichment models 

The core mass-radius relation is also relevant to self-enrichment 
models of old globular clusters, in which the cluster-forming core 
is often referred to as a 'protoglobular cloud'. Blue populations 
of globular clusters in elliptical galaxies show a 'blue-tilt', i.e. 
brighter clusters are redder than their fainter counterparts. This 
colour-magnitude relation is often interpreted as the imprint of the 
higher efficiency of more massive clusters to retain type II super- 
nova ejecta and to achieve greater metallicity. It must be borne 
in mind, however, that such a conclusion sensitivel y depends on 
the sl ope of the core mass-radius relation. While iMieske et al.1 
i nfer a positive slope for their predicted mass-metallicity 
relation, IParmentier & Gilmore ] 42001b predict that less massive 
clusters are more metal-rich. This discrepan cy arises because of 
different assumed core mass-radius relations. IMieske et al.l d2010l) 
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Table 2. Bound fractions F), om ^ integrated over a core mass function of slope — j8, upper bound m up and lower bound m; OM , = lOOMfi) , for the compact models. 
Other model parameters are as in Fig. [5] 
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build on either constant protoglobular cloud radii or constant pro- 
to globular cloud volume densiti es. On the other hand, the model 
of IParmentier & Gilmore I <200lh . applied to the Galactic Old Halo 
globular cluster system and developed before the 'blue-tilt' in el- 
lipticals was discovered, builds on pressure-bounded isothermal 
spheres for which 

t^core x t'core- Their conclusion that less mas- 
sive clusters are more metal-rich arises because type II supernova 
ejecta mix with a lower amoun t of primordial gas. Contrast i ng the 
models of iMieske etail J20 id) and lParmentier & Gilmorel feOOll) 
therefore illustrates that the slope of the globular cluster mass- 
metallicity relation predicted by self-enrichment models is a sen- 
sitive function of the assumed protoglobular cloud mass-radius re- 
lation. 



5 CONCLUSIONS 

We have investigated how the cluster-forming core mass-radius re- 
lation influences the resilience of clusters to an external tidal field 
during their violent relaxation, namely, while they dynamically re- 
spond to the expulsion of their residual star-forming gas. Owing 
to their high volume densities (nu2 — 10 4 ~ 10 5 cm -3 ), observed 
Galactic cluster-forming cores are virtually 'immune' to the exter- 
nal tidal field arising from an isothermal galactic halo with a cir- 
cular velocity V c ~ 220 km. s~ l . Following gas expulsion, however, 
clusters experience a severe density decrease through spatial ex- 
pansion and star-loss (cluster infant weight-loss). This renders ex- 
panded gas-free clusters more vulnerable to an external tidal field, 
possibly enhancing their infant weight-loss and likelihood of dis- 
solution compared to clusters evolving in a tidal-field-free environ- 
ment. 

To assess the impact of the core mass-radius relation 
on this issue, we have bu ilt on the /V-body model grid of 
iBaumgardt & Kroupa which provides the bound fraction 

Fbound °f stars m clusters at the end of violent relaxation as a func- 
tion of core SFE, gas expulsion time-scale TcExp I Across (in units 
of a core crossing-time T C ross) and tidal field impact r n /r t , where 
r/j and r, are the half-mass radius and tidal radius of the embed- 
ded cluster, respectively. Assuming a given external tidal field, SFE 
and gas expulsion time-scale, we have estimated the bound frac- 
tion Fbound as a function of the mass of cluster progenitor cores for 
various core mass-radius relations r core = X m core- We have consid- 
ered the cases of constant core surface density (5 = 1 /2), constant 
core volume density (S = 1/3) and constant core radius (8 = 0), 
which we refer to the S core , p core and r core models. For each slope 
8, we have tested two different normalizations %, which we re- 
fer to the 'compact' and 'loose' models (see TablefTJ. Loose mod- 
els describe molecular cores mapped in C 18 emission line. These 



cores are sometimes void of YSOs. In contrast, compact models 
describe dense molecular cores selected for their star-formation ac- 
tivity. The compact core mass-radius relations therefore constitute 
a better proxy to cluster-formation conditions than their loose coun- 
terparts. 

We have shown that constant surface density cores and con- 
stant radius cores result in the preferential removal of high- and 
low-mass clusters, respectively. In contrast, constant volume den- 
sity cores robustly lead to mass-independent cluster infant weight- 
loss (Fig. O. This is because constant surface density cores are 
characterised by a decreasing volume density - hence growing vul- 
nerability to an external tidal field - with increasing core mass, 
while the oppposite is true for constant radius cores (Fig.[4]l. 

Building on the compact mass-radius relations (Y. core — 
0.5 g.cm~ 2 , hhi = 6. 10 4 cm~ 3 and r core = 0.3 pc), these results are 
further quantified. At a galactocentric distance D„„; ~ 4 kpc in an 
isothermal galactic halo with a circular velocity V c = 220 km. s^ 1 , 
constant surface density cores preclude the formation of bound 
clusters more massive than a few 1O 4 M0 as a result of the greater 
vulnerability of more massive cores to the external tidal field. In 
contrast, the tidal field impact for constant volume density cores is 
independent of their mass and weak, that is, these cores respond to 
gas expulsion as if evolving in a tidal-field-free environment (see 
Figs|2][3]and|4j- The presence of massive star clusters (> 10 6 Mq) 
in the old haloes of elliptical and spiral galaxies, as well as in 
starbursts and galaxy mergers, may therefore constitute a hint that 
cluster-progenitor cores are not characterised by a constant surface 
density. Besides, almost all of the cluster-forming cores in the mid- 
dle and bottom panels of Fig.fJJhave ;i//2 > 10 4 cm -3 . Very massive 
molecular cores whose mean surface density is H C()re — 0.5 g.cm 
are not that dense. For instance, the mean number density of a 
core of mass 10 7 Mq and surface density L core = 0.5 g.cmT 2 is 
rifj2 — 750cm -3 only. Therefore, such cores may not result in the 
formation of massive clusters because of inefficient star formation 
in the first place. 

The shape of the mass function of young clusters is reported 
to be invariant through violent relaxation (albeit with a declin- 
ing amplitude), which points to mass-independent cluster infant 
weight-loss. In that respect, our analysis supports the hypothesis 
that cluster-forming cores have a constant volume density. When 
the SFE and gas expulsion time-scale are constant, cores of con- 
stant surface density produce mass-independent infant weight-loss 
for specific conditions only: the core mass range upper limit must 
be such that r/,/r, < 0.05. This can be obtained either through a 
weak tidal field (i.e. large r,) and/or a high surface density (to de- 
crease r/,), and/or a small upper limit on the core mass range as for 
observed molecular cores in our Galaxy (Fig.fTJ. 



Cluster-forming core mass-radius relation 13 



Figure [4] constitutes an efficient and straightforward tool to 
assess whether molecular cores produce clusters affected by an ex- 
ternal tidal field or not. Using Eq.Q/O] which provides the tidal field 
impact rf,/r t as a function of the core number density njj2 an d of 
the strength of an external tidal field (i.e. a galactic halo circular 
velocity V c and a galactocentric distance D ?a /), one can compare 
a given core mass-radius relation to lines of constant r/,/r, . As 
long as the core mass-radius relation stands below r/Jty = 0.05, 
the tidal field impact is weak to non-existent. As 0,/r, increases, 
cluster infant weight-loss becomes larger (see also Figs [2] and 
[3J. When r h /r, = 0.15-0.20, clusters survive only if the SFE is 
high and the gas expulsion time-sca le slow (e.g. SFE ^ 0.50 an d 
^GExp I Across ^ 3 when r/,/V/ — 0.20; lBaumgardt & Kroupa |2007| ). 
This is because high SFE and slow gas expulsion dampen the spa- 
tial expansion of an exposed cluster, thereby increasing its volume 
density and resilience to an external tidal field. 

We have emphasized that observationally inferred cluster- 
forming core mass-radius diagrams must be handled with caution 
as these can be imprints of the molecular tracer used to map molec- 
ular cores (e.g. C 18 vs. CS; see Fig. [TJ, or of the method used to 
measure the core radius and the mass enclosed within that radius 
(e.g. radius at the half-power beam width vs. radius w here the core 
densi ty profile has a characteristic value; see data of iMiiller et al.l 
2002, in middle and bottom panels of Fig. [I). 

Apart from the cluster mass function evolution, how the mass- 
radius relation of cluster-forming cores couples to an external tidal 
field to produce a tidal field impact upon expanded gas-free clusters 
also affects the star formation history of galaxies as reconstructed 
from the age distribution of th eir surviving star clusters (see also 
iMaschberger & Kroupa 1 120071) . Because constant surface density 
cores lead to the preferential removal of the most massive clusters, 
the bound fraction of stars integrated over the whole core mass 
function {Fj mmi ]) sensitively depends on adopted core mass func- 
tion parameters when 5 = 1/2 (spectral index — jS and upper mass 
limit m up \ see Table [2}. 

Independently of the tidal field impact, we note that the 
cluster-forming core mass-radius relation is also relevant to the du- 
ration of cluster violent relaxation and to the slope of the mass- 
metallicity relation of old globular clusters as predicted by self- 
enrichment models (see Section[4]l. 
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